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ABSTRACT 

The Spitzer Space Telescope has advanced debris disk science tremendously with a wealth of in- 
formation on debris disks around nearby A, F, G, K and M stars at 24 and 70 ^m with the MIPS 
photometer and at 8-34 fjtm with IRS. Here we present 160 /xm observations of a small sub-set of these 
stars. At this wavelength, the stellar photospheric emission is negligible and any detected emission 
corresponds to cold dust in extended Kuiper belt analogs. However, the Spitzer 160 /im observations 
are limited in sensitivity by the large beam size which results in significant "noise" due to cirrus and 
extragalactic confusion. In addition, the 160 ^m measurements suffer from the added complication 
of a light leak next to the star's position whose flux is proportional to the near-infrared flux of the 
star. We are able to remove the contamination from the leak and report 160 /im measurements or 
upper limits for 24 stars. Three stars (HD 10647, HD 207129, and HD 115617) have excesses at 160 
/im that we use to constrain the properties of the debris disks around them. A more detailed model 
of the spectral energy distribution of HD 10647 reveals that the 70 and 160 /im emission could be due 
to small water ice particles at a distance of 100 AU consistent with Hubble Space Telescope optical 
imaging of circumstellar material in the system. 
Subject headings: 



1. INTRODUCTION 

Surveys of nearby main-sequence stars with the Spitzer 
Space Telescope have produced a new sample of debris 
disks which will be studied for years to come. Since 
Foynting-Robinson (F-R) drag, radiation pressure and 
wind drag^ remove dust from the disk, the dust in debris 
disk systems must be replenished, presumably through 
collisions between larger bodies, possibly remnant plan- 
etesimals that were the building blocks of full-fledged 
planets. At 70 /im Spitzer MIFS photometry is sensitive 
to dust emission which is more than 5-10 times more 
luminous than that seen in our own solar system with 
excesses analogous to emission from our Kuiper belt oc- 
curring in ^--^15% of main sequence stars (Bryden et al. 
2006). At 8-14 /im Spitzer is sensitive to dust emission 
which is more than 1,000 times more luminous than that 
seen in our own solar system and excesses analogous to 
emission from our own asteroid the Kuiper belt are found 
to occur in roughly 1-2% of main sequence stars (Beich- 
man et al. 2006). This paper focuses on the detection of 
excesses at 160 /im thus probing the properties of Kuiper 
belt analogs around these stars. 

Most if not all, published Spitzer debris disk surveys 
have focused on either IRS spectra (14-34 /im) or MIFS 
24 and 70 /im photometry with several surveys reporting 
160 /im data for young stars (RebuU et al. 2008; Cieza et 
al. 2008). While the 24 and 70 /im results suggest that 
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160 /im observations should often show excess emission, 
Spitzer is much less sensitive at this wavelength due to 
confusion from extra-galactic sources in its large beam 
and to an artifact in the data known as "the leak" . This 
paper presents the results of our analysis of the 160 /im 
data taken as part of a Spitzer survey of nearby main 
sequence stars. Basic debris properties are derived from 
models of the spectral energy distribution of those stars 
with excesses at 160 /im. Section 2 briefly discusses the 
target sample considered in the study, § 3 reviews the 
reduction of the images, § 4 discusses the analysis used to 
remove the leak and determine the 160 /im photometry, 
§ 5 summarizes the statistics of our findings and SED 
models of the data and § 6 reviews the implications of 
our results on the overall properties of the disks around 
mature stars. 

2. SAMPLE 

The primary sample of stars analyzed in this paper 
originates from the collection of mature FGK stars ob- 
served with MIFS at 24 and 70 /im (Beichman et al. 
2005; Bryden et al. 2006, Program Identification (FID) 
#41). To aid with the analysis of the leak, we also in- 
clude 13 bright stars observed at 160 /im as part of an 
M dwarf program (G. Rieke, FID 52) or in a MIPS cal- 
ibration program (B. Latter, FID 1709). There were no 
criteria on the spectral type of the stars used to model 
the leak since they are meant to be bright FSF templates 
and the 160 /im flux originates from the Rayleigh Jeans 
portion of their blackbody spectra. Tables [1] and [2] list 
the science targets and calibration stars. 

3. OBSERVATIONS AND DATA REDUCTION 

All Spitzer observations were made between December 
2003 and March 2006 with integration times of 3 and 10 
seconds. The data reduction is based on the DAT soft- 
ware developed by the MIPS instrument team (Gordon et 
al. 2005) and is similar to that performed on the 24 and 
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70 /im data (Beichman et al. 2006; Trilling et al. 2007; 
Bryden et al. 2006). We used images processed beyond 
the standard DAT software and then mosaicked from in- 
dividual frames with half pixel subsampling (Stansberry 
et al. 2007). The resulting pixel scale is 8" pixel" ^. 
The calibration factor of 0.968 /iJy arcsec"^ DN~^ s~^ 
and the appropriate aperture corrections for the 160 fim 
MIPS images were taken from Stansberry et al. (2007). 

4. ANALYSIS 

The 160 fim data suffers from an imaging artifact 
(a.k.a. the "leak") that is caused by the reflection within 
the instrument of near-infrared light onto the MIPS de- 
tector (Stansberry et al. 2007, see Figure [T]). The leak 
becomes an issue when looking for excesses around bright 
stars since the near-IR photosphere is bright compared 
to the weak 160 /j,m signal. The leak lies ~4 pixels to 
the +X side of the star's position on the detector and is 
close enough to hinder accurate 160 /xm flux determina- 
tions with aperture photometry. 

Since the leak is due to the near-infrared light from the 
star, we can predict the level of the leak emission based 
on the known near-IR brightness of the star. We utilized 
archived Spitzer observations of a sample of stars which 
have no 160 /im or shorter wavelength excess and which 
are bright in the near-infrared (see Table [2]). These stars 
vary in brightness from Ks = —1.39 (Sirius) to Kg = 4.23 
(HD 166620). We assume that the dominant source of 
flux in the 160 /im image near the position of the star is 
the leak. The offset of the 160 /im leak from the position 
of the star, determined within a fraction of an arcsecond 
from the contemporaneous 24 /im MIPS image, has a po- 
sition and dispersion of 4.5±1.5 pixels in the X direction 
and -O.libl.O pixels in the Y direction (see Figure [2|). 
The spread in the leak offset is presumably due to differ- 
ences in the observing parameters of the targets which 
have been collected from a number of different AORs. 

To estimate the flux of the leak we used aperture pho- 
tometry with an aperture radius of 3 pixels (24") and a 
sky annulus of 3-5 pixels (24-30"). Figure [3] plots the 160 
fj,m leakage versus the 2MASS Ks flux for the leak-only 
stars in our sample. There is a clear dependence between 
the 160 ^m leakage and K band flux with a correlation 
coefficient of 0.98. The best fitting line through the data 
is Fi^ak[mjy] = (10.5 ± 1.0)Fj.2^m[jy] + 26 ± 15. There 
is no correlation between the leak flux or offset and the 
J-K color of the star. 

To remove the leak from our target sample, we create 
a "leak" Point Spread Function (PSF) by scaling and av- 
eraging the four brightest leak-only images (Sirius, Pro- 
cyon, HD 80007, and HD 27442). The derived "leak" 
PSF is shifted to the appropriate position, scaled to the 
predicted flux and subtracted from the image. Figure 2] 
(middle pane) shows the result of this subtraction for HD 
10647 and HD 207129. In some cases the leak was not 
fully removed using the scaled PSF. This could be due to 
large errors in the K5 band flux due to saturation in the 
2MASS images as suggested by the K5 error bars shown 
in Figured 

As an alternative method to remove the leak, we em- 
ploy a cleaning method which subtracts a scaled version 
of the leak iteratively in order to reduce the signal level in 
the leak region to match that in the surrounding back- 
ground. First, a circular region with a radius of four 



pixels at the position of the leak is declared the clean- 
ing area. During the cleaning process, the highest point 
within the cleaning region is identified and the leak PSF 
is scaled to 10% of the value of the pixel at that position. 
This scaled leak PSF is subtracted at this position and 
then the process is repeated until the RMS noise within 
the cleaning area is at the same level as that in an an- 
nulus around the position of the star. Figure [4] (bottom 
pane) shows the results of the iterative method for HD 
10647 and HD 207129. Figure [D shows the results of the 
leak removal process for HD 142860 which does not have 
significant excess emission at 160 ^m. These cleaned im- 
ages show that the leakage is getting reduced to the level 
of the background using the iterative method. We use 
the results from both the scaled leak and the iterative 
cleaning methods in determining the uncertainty in the 
160 /im fiux. 

We performed aperture photometry on the leak- 
removed and raw images using an aperture of 4 pixels 
(32") in radius and a sky anmflus of 7-8 pixels (56-64"). 
Figureinishows a histogram of the flux densities measured 
within the aperture before (empty) and after (diagonal) 
removal of the leak for all stars in the 160 /im sample. 
Removing the leak results in an 40% improvement in the 
flux sensitivity of the data based on the average values 
of these two distributions. For those stars with no signif- 
icant 160 /im excess, we estimate an upper limit to the 
160 /im photometry based on the flux that would be nec- 
essary for an excess to be detected given our detection 
criteria and background noise estimate for each image. 
Table 2] lists the upper limits for those stars with 70 /im 
excesses and no 160 /^m excesses. Most of these upper 
limits provide a more stringent constraint on the range of 
temperatures expected for the cold dust population that 
would be estimated from solely the 70 /im detection. 

5. RESULTS 

For this study, we use the criterion that a star has a 
160 /im excess if the flux at the position of the star after 
leak removal and subtraction of the stellar photosphere 
is 3 (7 above the background noise in the image. To de- 
termine the level of background noise, we flrst correct for 
the flux gradient apparent over the length of the MIPS 
160 array. The gradient is removed by fltting a line to 
pixels at both ends of the array and using the slope to 
create a gradient which is subtracted from the image. 
The background noise is then estimated from the stan- 
dard deviation of the counts in eight apertures placed 
14 pixels (122") to the right {+X) of the star. Using 
these criteria, only three stars have significant 160 /im 
emission after subtraction of the leak: HD 10647, HD 
207129, and HD 115617. Table [3] lists the 160 /im flux 
densities estimated with aperture photometry from the 
non-leak subtracted image, the image with the K-band 
scaled leak subtracted and the image with the leak re- 
moved iteratively. We also list the background noise flux 
density and the flux density of the stellar photosphere at 
160 /im. Unlike the 24 and 70 /im studies, all detections 
at 160 /im are dominated by the excess emission above 
the stellar photosphere. 

Figure[8]plots the minimum detectable disk luminosity, 
Ltj/L*, estimated from the 160 /im excess as a function 
of stellar effective temperature, T*. This assumes a sin- 
gle disk temperature which is determined by setting the 
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blackbody emission peak at 160 fj,m corresponding to a 
dust temperature of 23 K. The minimum disk luminosity 
is then the ratio of the total flux determined from this 
blackbody to the stellar photospheric flux (see Equation 
1, where Figo is the 160 fim flux density and Figg,* is the 
flux density of the stellar photosphere at 160 /im). Most 
of these luminosity limits (see Table [S]) are an improve- 
ment over those provided by 70 /im limits from previous 
studies (Bryden et al. 2006). 

Ldustmin „ f ^^^^^\^ -^160 n 

= 8 X 10 I — - — I (1) 



Using these L^/L* values, we can provide an estimate 
of the dust mass assuming 1 /im radius grains with a 
density of 3.3 g cm~'^ at a distance of 100 AU (Jura et 
al. 1995, see Table [5]). 

For those stars with 70 /im excesses and 160 /im upper 
limits, we can estimate a lower limit for the tempera- 
ture of the dust component responsible for the excess 
at these two wavelengths by fitting a blackbody curve to 
these two data points (see Table [5]). Unfortunately, these 
temperature limits are more indicative of the variations 
in the background levels of the 160 /im images than the 
properties of the dust. As an example. Figure [7] shows 
spectral energy distributions of HD 17925 and HD 82943 
with a blackbody fit to their 70 /im excess and 160 /im 
upper limit. Equation 2 shows that we can also estimate 
a limit on Ljj/L* (also in Table [6]) in the same manner as 
for the 160 /im detections but with the 3 a flux limits. 
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where T, and Figo,* are the stellar effective temperature 
and flux, respectively, and Nigo is the one a flux den- 
sity upper limit. Finally, we can estimate the smallest 
amount of dust that could be present to have no 160 /im 
detection. Using the L^/L* lower limits, we can provide 
a lower limit to the dust mass using the same assump- 
tions for those stars with 160 /im excesses. Table [5] lists 
these dust masses as fractions of an Earth mass. 

5.1. Spectral Energy Distributions 

Here we employ simple SED models to fit the 24, 70 
and 160 /im MIPS photometry and IRS spectra from 
previous studies of the three stars with 160 /im excesses. 
The 24 and 70 /im photometry comes from Trilling et al. 
(2008) and the IRS data come from either Beichman et 
al. (2006) or Lawler et al. (2009). AU of the stars with 
160 /im detections have disks which are resolved at 70 /im 
(HD 10647, 106 AU (radius); HD 207129, 144 AU; and 
HD 115617, 97 AU; Bryden et al. 2009 in prep). The 
70 /im sizes for HD 10647 and HD 207129 both agree 
with the Hubble Space Telescope (HST) measurements 
(Stapelfeldt et al. 2007) which detect a relatively nar- 
row toroidal structure at 90-140 AU using the Advanced 
Camera for Surveys (ACS). 

5.2. Silicates or Black Body Dust Models 

To model the debris disk we place the dust in an op- 
tically thin annulus around the star and utilize the for- 
malism of Su et al. (2005). The dust temperature varies 
with the distance, r, from the star as determined from 



radiative equilibrium. The free parameters in the SED 
model fits are the inner and outer radius of the dust 
annulus, and Tq, and the dust density at the inner 
radius which we later convert to L^/L* and then the 
total dust mass in the system. Also, we make assump- 
tions of the dust composition which affects the shape of 
the SED We choose to fit the Spitzer data with mod- 
els incorporating either blackbody grains, small silicate 
dust grains (a = 0.25 /im) or large silicate dust grains 
(a = 10 /im), with emission coefficients taken from Laor 
& Draine (1993). The smallest silicate grains used in our 
models are bigger than the blowout radius of the grain 
which depends on the luminosity of the star. Since 10 /im 
dust grains are primarily affected by Poynting-Robinson 
drag, we model this disk component with a constant sur- 
face density. The 0.25 /im grains, on the other hand, are 
primarily removed from the disk by radiation pressure 
resulting in a surface density with a r~^ fall off (Su et 
al. 2005). If we include the additional constraint that 
the 70 /im MIPS photometry must be modeled with a 
dust component with an inner radius corresponding to 
the resolved 70 /im radius, then the Spitzer data is best 
fitted with two separate dust components with different 
dust species. 

The properties of the dust derived from the best fits 
to the Spitzer data (the inner and outer radii of the dust 
annulus, L^/L*, and the dust composition) are listed in 
Table [5] along with derived dust properties such as dust 
temperature and dust mass. Figure [9] shows the best fit- 
ting SED models for all three stars with 160 /im excesses. 
The errors on the variables used to fit the SED, the in- 
ner and outer dust radius and inner dust density of each 
dust component, were estimated from the change in their 
value which produced a variation of one in the (Press 
et al. 2002). 

5.3. Icy Dust Models 

A more thorough investigation of the SED of HD 
10647, based on the primitive solar system dust models 
of Lisse et al. (2007, 2008) reveals that, after model- 
ing the IRS 5-35 /im warm dust signature, the 70 and 
160 /im fiux densities can be modeled by a component 
of '^30 K water ice particles at 100 AU from the star. 
An enhanced peak in the 70 /im passband occurs due to 
emission from a strong water ice vibrational feature at 
65 /im (see Figure [TU)) . While additional contributions 
to the dust emission in the IRS portion of the spectrum 
could be due to silicates, our data set cannot place statis- 
tically significant constraints on the abundance of such 
materials. These SED models do, however, showcase the 
improvement in our ability to learn about the dust with 
continuous spectral coverage and direct spatial informa- 
tion. A distance of ~100 AU from the primary is con- 
sistent with the 90 - 140 AU location of dust detected in 
the optical with HST/ ACS (Stapelfeldt et al. 2007). 

In addition to the Spitzer observations presented here, 
thermal emission from the HD 10647 circumstellar dust 
has been studied most recently at 870 /im using the sub- 
mm LABOCA array on the 12-m APEX telescope in 
Chile (Liseau et al. 2008). In order to fit the 870 /im 
flux point of 39±8 mJy for this star, the dominant (by 
surface area) water ice particle sizes must be on the or- 
der of 1 mm. However, as stated in Liseau et al. (2008), 
the morphology of this cold component is very broad, on 
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the order of 600 AU FWHM. It is most likely that the 
material emitting the bulk of the 870 /xm emission is at 
17 K, and is unassociated with the material constrained 
to the 100 AU torus. We are unable to fit the 870 /im 
data point with the water ice SED component at 30 K 
which is consistent with the hypothesis of an additional 
17 K dust component at a distance of ^ 300 AU. 

6. DISCUSSION 

For this sample we only consider those stars in the orig- 
inal FGK survey (PID41, Beichman et al. 2005; Trilling 
et al. 2007) with 160 fj,m data. While all stars with 
160 fim excesses also have 70 /im excesses, not all stars 
with 70 /im excesses have 160 fj,m excesses. Out of 35 
stars in the FGK sample, two have 160 /im excesses. 
The other one, HD 10647, was included in a separate 
survey with different goals. Twelve of the 160 /zm non- 
detections have 70 fim excesses. One (HD 10647) of the 
five stars with known RV planets has a 160 /im excess 
(20%), compared to 2 out of 31 (6%) of the stars with 
no known planets. This sample is smaller than the 24/70 
/Ltm samples so these percentages suffer from small num- 
ber statistics. Given the 160 /um sensitivity limits due to 
confusion noise and the leak, we are generally not sen- 
sitive to dust much colder than 25 K and 0.1x10"^ M0 
(see Table E]). 

A collection of icy dust with a mass of '^4x10"'* 
at 30 K around HD 10647 is very reasonable with respect 
to known structures in the solar system. Dust created by 
interactions in the Saturnian and Uranian systems (such 
as formed the system's rings) will equilibrate at 60-90 
K. Dust created by coUisional fragmentation of bodies 
in the Kuiper Belt (Davis & Farinella 1997; Stern 1995) 
produces originally toroidal morphologies at local tem- 
peratures of 35-45 K. The location of the dust is far re- 
moved from the orbit of the reported planet HD 10647b 
with semi-major axis of 2.1 AU (Mayor et al. 2003) and 



is not likely to be dynamically coupled to it. 

7. CONCLUSIONS 

We have analyzed the 160 /im subset of a survey of 
nearby main sequence stars with the Spitzer IRS and 
MIPS instruments. After careful analysis of the 160 fim 
images which suffer from an anomalous flux contaminant, 
we have found three stars (HD 10647, HD 207129 and HD 
115617) with significant excesses in this passband. All of 
these stars also have excesses at 70 /im and the addi- 
tion of the 160 /im photometry allows for a constrained 
estimate of the temperature of the dust as well as its 
composition. The lack of 160 /im emission for the re- 
maining stars also sets an improved limit to the amount 
of cold dust for those sources with 70 /im excesses. In 
the end, the sensitivity floor of the 160 /im images is set 
by the background noise due to extragalactic sources. 

The Hershel Space Telescope, scheduled to launch in 
2009, will have about four times better resolution than 
Spitzer and will be '-^lO times more sensitive at compa- 
rable wavelengths suggesting we might be able to resolve 
those disks detected at 160 /im. Having the intrinsic size 
of the debris disks as well as robust fluxes at multiple 
wavelengths would further constrain the SED models. 
Those stars with 70 /im excesses and no 160 /im detec- 
tions are ideal targets as well since there is dust emission 
at longer wavelengths not detectable with Spitzer. 

The research described in this publication was carried 
out at the Jet Propulsion Laboratory, California Insti- 
tute of Technology, under a contract with the National 
Aeronautics and Space Administration. This publica- 
tion makes use of data products from the Two Micron 
All Sky Survey, which is a joint project of the University 
of Massachusetts and the Infrared Processing and Analy- 
sis Center /California Institute of Technology, funded by 
the National Aeronautics and Space Administration and 
the National Science Foundation. 
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TABLE 1 
Stars in the 160 ;um Sample 



Star 


SpTy 


Ks 


Distance [pc] 


PID" 


tiL) ±DO 


ivu ve 


O.io 


1 Q '7n 


A 1 
41 


HD 10478 


KIV 


5.20 


7.47 


41 


TT"p\ 1 ridAT b 

rlJJ lUo4/ 


r y V 


5.52 


17.35 


716 


HD 13445 


KIV 


6.17 


10.91 


41 


HD 17925 


KIV 


6.00 


10.38 


41 


HD 33262 


F7V 


4.72 


11.65 


41 


HD 33636 


GO 


7.06 


28.69 


41 


tiu 37cSy4 


KIV 


6.23 


12.24 


41 


HD 52265*' 


GOIII-IV 


6.30 


28.07 


41 


HD 72905 


G1.5Vb 


5.65 


14.27 


41 


HD 76151 


G2V 


6.00 


17.09 


41 


HD 82943*" 


GO 


6.54 


27.46 


41 


HD 114783*' 


KO 


7.57 


20.43 


41 


HD 115617 


G5V 


4.74 


8.53 


41 


HD 117176 


G5V 


5.00 


18.11 


41 


HD 118972 


K2V 


6.93 


15.61 


41 


HD 128311 


KOV 


7.51 


16.57 


41 


HD 145675 


KOV 


6.67 


18.15 


41 


HD 149661 


K2V 


5.76 


9.78 


41 


HD 177830*' 


KO 


7.18 


59.03 


41 


HD 185144 


KOV 


4.70 


5.77 


41 


HD 190007 


K4V 


7.48 


13.11 


41 


HD 207129 


GOV 


5.58 


15.64 


41 


HD 219134 


K3V 


5.56 


6.53 


41 


HD 220182 


KIV 


7.36 


21.92 


41 



PID = program identification number'' Stars with known planet(s). 



TABLE 2 
Stars I'sed foh tiii; Leak PSF 



Star 


SpTy 


V 


Ks 


Distance [pc] 


PID" 


HD 3651 


KOV 


5.80 


4.00±0.04 


11.11 


41 


HD 4628 


K2V 


5.75 


3.68±0.27 


7.46 


41 


HD 27442 


K2IVa 


4.44 


1.75±0.22 


18.23 


41 


HD 80007 


A2IV 


1.70 


1.49 ±0.24 


34.08 


1709 


HD 88230 


K5V 


6.61 


2.96±0.29 


4.87 


41 


HD 142860 


F6IV 


3.85 


2.70±0.31 


11.12 


41 


HD 166620 


K2V 


6.37 


4.23±0.02 


11.10 


41 


HD 191408 


K3V 


5.31 


3.01±0.60 


6.05 


41 


HD 203608 


F6V 


4.22 


2.97±0.25 


9.22 


41 


HD 209100 


K4.5V 


4.69 


2.24±0.24 


3.63 


41 


HD 216803 


K4V 


6.48 


3.81±0.24 


7.64 


41 


Procyon 


F5IV-V 


0.34 


-0.66±0.32 


3.50 


52 


Sirius 


AlV 


-1.47 


-1.39±0.21 


2.64 


52 



PID = program identification number 



TABLE 3 

Flux Densities for Stars with 160 nvi excesses 



Star 


Raw"" 
mjy 


Excess - scaled leak method*' 
mJy 


Excess - iterative method*' 
mJy 


Noise"^ 
mJy 


Significance"* 
mJy 


Stellar Photosphere*^ 
mJy 


HD 10647 


498 


462 


451 


50 


9.2 


3.17 


HD 207129 


185 


152 


158 


20 


7.9 


3.45 


HD 115617 


318 


141 


89 


20 


4.5 


10.56 



^ The raw flux represents the combination of the flux from the leak, the stellar photosphere and excess thermal emission.'' The "scaled leak 
method" and "iterative method" represent two different forms of analysis utilized to removed the leak flux prior to performing aperture 
photometry at the position of the star. The difference in the two methods are described in Section 4. The background noise is estimated 
from the standard deviation on the flux from eight apertures placed 35 pixels in the +x direction away from the star. The significance 
given here is defined as the excess thermal flux (sky subtracted) from the iteratively cleaned image over the background noise." Value of 
the stellar photosphere at 160 /xm. This has been subtracted from all of the excess fiux values. 
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TABLE 4 
Stars With No 160 fj,u Excess 



Star 


Flux Density Limit [mjy] 


70 /im excess? 


HD 166 


<30 


yes 


xlU iryzo 


< /O 


yes 


JlU OoZOZ 




yes 


rlJJ ooDoD 


<oz 


yes 


TT T\ O'JQCiA 


<7o 


yes 


TT cooi^c; 
rlJJ OZZOO 


<523 


yes 


TT vonnK 

rijj i zyuo 


<55 


yes 


TTr~\ Vi^i cri 
rlJJ iOiOi 


<iou 


yes 


TTT~\ QOn/IQ 


<iOU 


yes 


upl 1171 


<,.oUU 


yes 


HD 118972 


<110 


yes 


HD 128311 


<83 


yes 


HD 10476 


<90 


no 


HD 13445 


<20 


no 


HD 114783 


<60 


no 


HD 145675 


<47 


no 


HD 149661 


<87 


no 


HD 177830 


<343 


no 


HD 185144 


<372 


no 


HD 219134 


<1334 


no 


HD 220182 


<60 


no 



TABLE 5 
Dust Model Results 



Star 






^ dust 


Ld/L* 


Dust Mass" 


Dust Composition'' 




AU 


AU 


K 


10-5 


10-4 M® 




HD 10647 


6±0.5 


150±150 


70-40 


21.1 


3.43 


blackbody 




97±10 


140±50 


33-27 






silicate - 0.25 /im or water ice 


HD 115617 


96±5 


195±100 


55-45 


2.6 


0.50 


blackbody 




120±20 


220±10 


24-19 






silicate - 0.25 /im 


HD 207129 


35±1 


45±10 


48-44 


6.6 


0.16 


silicate - 10 /im 




144±1 


200±10 


24-12 






blackbody 



^ The dust mass is estimated from the L^j/L* ratio, a dust density of 3.3 g cm , a dust radius of 1 /im, and a disk radius of 100 AU as 
defined in Jura et al. (1995).*^ Here, we list the assumed composition and radius (when necessary) of the grains used in the models. We 
assume either black body grains or silicate grains. The silicate grains can have a radius of 0.25 jxra. (small grains) or 10 jxra. (large grains). 
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TABLE 6 

Dust Physical Properties for Stars with 160 ixu Upper Limits 



Star 


Dust Temperature" 




Dust mass"^ 




K 


10-5 


10-" Me 


HD 166 


>104 


<0.16 


<0.33 


HD 17925 


>42 


<0.33 


<0.68 


HD 33262 


>27 


<0.20 


<0.41 


HD 33636 


>34 


<0.73 


<0.22 


HD 37394 


>23 


<0.80 


<1.51 


HD 52265 


>21 


<2.90 


<6.02 


HD 69830 


>21 


<1.83 


<0.90 


HD 72905 


>36 


<0.17 


<0.35 


HD 76151 


>25 


<0.64 


<1.33 


HD 82943 


>43 


<1.11 


<2.30 


HD 117176 


>24 


<0.50 


<1.03 


HD 118972 


>29 


<1.04 


<2.15 


HD 128311 


>27 


<1.28 


<2.66 


HD 10476 




<0.19 


<0.39 


HD 114783 




<1.03 


<2.13 


HD 13445 




<0.09 


<0.19 


HD 145675 




<0.33 


<0.69 


HD 149661 




<0.29 


<0.60 


HD 177830 




<31.8 


<65.9 


HD 185144 




<0.46 


<0.96 


HD 219134 




<3.51 


<7.28 


HD 220182 




<0.87 


<1.S0 



* The dust temperature listed here represents the blackbody temperature that reproduces the ratio of the 70 /xm flux and 160 //m 
limit.'' These values of L^j/L* come from Equation 2.*^ Upper limits on the dust mass assume the L^/L* ratio from the 160 /im limit, a 
dust density of 3.3 g cm-*, a dust radius of 1 /im, and a disk radius of 100 AU (Jura et al. 1995). 



HD10647 




Arcseconds 

image of HD 10647 at 160 fim showing the position of the star relative to the leak. 
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Fig. 2. — Plot of the variation in the offset in pixels between the leak ccntroid and the position of the star as estimated from the star's 
24 /im position. 
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10 100 1000 10000 

Kflux [Jy] 

Fig. 3. — Plot of the K band flux from 2MASS versus the flux of the leak estimated from a sample of stars with 160 fj,m flux below Spitzer 
sensitivities. The plot shows a clear correlation between near-infrared flux and peak flux allowing for the potential of removing the leak flux 
using this correlation as a guide. The best fitting line which is flt through all data points is Fisaklmjy] = (10.5 ± l-0)i^2.2/im[ji/] + 26 ± 15. 
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Fig. 5. — MIPS 160 fim images of HD 142860 before (top) and after the leak was removed using the scaled leak (middle) and iterative 
(bottom) cleaning processes. In this case, there is no detection at 160 /im but its clear that the leak has been fully removed from the image 
using the iterative method. The plate scale of the images is 8" pixel^^. 
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Fig. 6. — Histogram of the background flux density at the positions of the 160 fim stars before (empty) and after (diagonal) leak removal. 
Removing the leak results in an 40% improvement in the flux sensitivity of the data based on the average values of these two distributions. 




Fig. 7. — SEDs of HD 17925 and HD 82942 which do not have significant 160 fim excesses. The upper Hmit in 160 fim flux density can 
still be used to constrain the properties of the disk. 
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Fig. 8. — Plot of the limits on L^j/L, for all stars in the 160 /im sample as a function of stellar effective temperature. The crosses, +, 
those stars with detected 160 fim flux. Those stars with non-detections have upper limits in their L^/L, values provided by Equation 
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Fig. 9. — SEDs of those stars with 160 fim excesses. In all cases, these stars also have excesses at 70 microns allowing for a fit to the 
temperature of the dust. The dust is placed in an annulus around the star with a number density fall off dependent on the dust removal 
mechanism. The additional constraint of the extent of the dust at 70 ^m and the IRS spectra, results in the best fit (solid line) being 
achieved with two separate dust annuli (dotted and dashed lines) of varying composition. The dust temperatures (20-70 K) were estimated 
assuming radiative equilibrium and depended on the dust composition, stellar luminosity and distance of the dust from the star). The 
properties of the dust derived from these SED model fits are listed in Table [5] The triangles are IRAS (25, 60 fim) fluxes and the squares 
are the MIPS data. 
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Fig. 10. — Plot of an additional model of the Spitzer data for HD 10647. This model uses 30 K water ice particles to model the 70 and 
160 /im emission component (dotted) and ~70 K silicates to model the IRS data (dashed). The model has been normalized to the MIPS 
24 /im flux density. 



